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Abstract 

We reformulate MHV scattering amplitudes in 4D gauge theory and supergravity as 
correlation functions of bilinear operators in a supersymmetric gaussian matrix model. The 
model retains the symmetries of an S'^ of radius i and the matrix variables are represented 
as linear operators acting on a finite-dimensional Hilbert space. Bilinear fields of the model 
generate a current algebra. In the large double scaling limit where ipi ~ i/\/N is held 
fixed, there is an emergent flat 4D space-time with a built in short distance cutoff. 
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1 Introduction 



Recent progress in the study of scattering amplitudes in supersymmetric gauge theory and 
gravity has revealed surprising structures that hint at the existence of deep new principles 
[THTT]. This development was triggered by, and has gone hand in hand with, the exploration 
of the various dualities between gauge theory, gravity and string theory - most notably the 
AdS / CFT correspondence, and the proposed reformulation of A/" = 4 SYM amplitudes via 
twistor string theory [12], P^fT^ (see [13] for an early review). 

While influential in the beginning, the twistor string program has been relatively dor- 
mant in recent years, in large part because of the realization that the theory contains 
conformal supergravity, and by the recognized inconsistency of the latter. In spite of this 
apparent roadblock, the twistor string is clearly a natural and beautiful idea that deserves 
to find its place among the realm of consistent string theories. Rather than a fatal weak- 
ness, the unexpected emergence of gravity could count as an extra motive for trying to get 
to the bottom of its true significance. 

In the accompanying paper [TH], we propose a new interpretation of twistor string theory 
that potentially avoids the trap that led to its apparent inconsistency. The new idea 
is to view the strings as emergent low energy degrees of freedom of a pure holomorphic 
U{NNc) Chern-Simons gauge theory in the presence of a U{N) background flux with a 
large instanton number The twistor strings then arise as collective modes of the 

instantons. Their dynamics is captured by an effective theory that contains a U{Nc) gauge 
field A, also with a hCS action, living on a non-commutative twistor space. In addition, 
the instantons give rise to a pair of defect modes Q and Q. As shown in [TB], the low energy 
gauge field and defect modes naturally combine into a matrix model that encompasses the 
ADHM construction of instantons. The matrix model is initially formulated at finite 
and comes with a length scale i given by the radius of an S*^. See [23] for earlier work on 
potential connections between matrix models and twistor string theory. 

As a concrete output of this new proposal, we have identified a simple gaussian large A^ 
matrix model given by an action of the form: 



where Q and Q are finite matrices, and is a non-commutative covariant derivative, that 

^ Our choice to consider the large instanton limit was stimulated by previous proposals that 
formulate theories of quantum gravity in terms of large N matrices as for example in |18H22j . 
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is, a linear operator acting on Q. This gaussian model is obtained by starting from the 
interacting ADHM matrix model introduced in [16], and working with respect to a fixed 
background for the non-commutative gauge field A. Our goal in this paper is to establish 
that this truncated model describes the MHV sector of an emerging space-time theory. 

The gaussian matrix model (II. ip enjoys a number of remarkable properties. The model 
comes with an 50(5) symmetry, that refiects its relation with an underlying 5*^ space- 
time geometry. Even at finite A^, the matrices acquire a natural geometric interpretation 
as holomorphic functions on twistor space. As a result, the basic link (see e.g. [2^125]) 
between twistors and space-time physics is preserved. The fiat space continuum limit is 
obtained by taking a double scaling limit N ^ oo with: 




held fixed. In this limit, we will identify spin one excitations for a u{N^ gauge theory 
as well as spin two excitations corresponding to deformations of the space-time geometry. 
This motivates an identification of the short distance cutoff Ipi with the Planck length. 

The basic link we establish in this paper is that in this double scaling limit, states with 
specified momentum and helicity are represented by currents J made from bilinears in the 
matrix variables. Correlation functions of multiple Ji insertions compute amplitudes of the 
4D theory via the correspondence: 



This link between currents and space-time fields is somewhat similar to the AdS/CFT 
dictionary, and to the map between target space fields and vertex operator of the string 
worldsheet theory. MHV gauge theory amplitudes are computed via a u{N^ current al- 
gebra, and naturally reproduce the Parke- Taylor formula [2S]- Moreover, we find that the 
matrix model does not give rise to conformal gravity amplitudes, but rather, in the strict 
large A^, Ipi — > limit, produces an effective action that matches the generating functional 
of MHV amplitudes in A/" = 4 SYM theory introduced by PTIEH]. 

A striking feature of the matrix model is that it also allow for other geometric currents, 
that, as long as Ipi is kept finite, generate complex structure deformations of the non- 
commutative twistor geometry. The properties of these currents and deformations are 
strongly reminiscent of the non-linear graviton construction of Penrose [21]. Motivated by 
this relationship, we will compute correlation functions of these geometric currents and 
quite remarkably, we will find that they reproduce the BGK amplitude for MHV graviton 




(1.3) 
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scattering pO]. In this sense, we can identify a natural algebra associated with MHV gravity 
amplitudes, which is also reminiscent of BCJ [31] (see also [32]). The final representation 
(and several geometric and combinatorial steps in the derivation) of the amplitude takes 
the same form as the formula derived in 



{nl)^ 1 '}-^[k\pi + ... +pk-i\n) 



Jn-l)(n-ln)(nl)C(n)Al {kn) ^ ^(2,...,n-i)- 

(1.4) 



/ [ln — \){n — ln){nl)C{n)^^^ 

where the sum P2,...,n-i is over all permutations of the plus helicity gravitons. Here, n = 
\^167tGn, and C{n) = (12) (23) • ■ ■ {n — ln){nl) is the usual Parke- Taylor denominator. 

Our plan in this paper will be to start from the gaussian matrix model, and to show 
how features of the 4D space-time theory are built up. This can be viewed as a "bottom 
up" perspective on the matrix model and its connection to a potential theory of emergent 
space-time and gravity. In the companion paper [16], we provide a more UV motivated 
perspective on the proposal. 

The rest of this paper is organized as follows. In section[2]we introduce the basic features 
of the gaussian matrix model which follow from [16l[T7] . Section [3] studies the symmetries 
and currents of the matrix model. In section H] we provide a space-time interpretation for 
the matrix model. This will motivate the conjecture that correlators of the matrix model 
are connected with 4D physics. In section [5] we study correlators for a fuzzy CP^. This is 
of interest in its own right, and will serve as a springboard for the full computation of the 
matrix model correlation functions. We formulate the calculation of scattering amplitudes 
in section [6], and show that MHV gluon correlators are naturally reproduced from a u{Nc) 
current algebra. In section [7] we compute MHV graviton scattering amplitudes. Section [8] 
contains our conclusions. Some additional technical details are included in the Appendices. 



2 Gaussian Matrix Model 



In this section we introduce the supersymmetric gaussian matrix model. The matrix vari- 
ables are represented as linear operators acting on a finite dimensional Hilbert space, ob- 
tained by quantizing A/" = 4 supersymmetric twistor space CP'^'^. 
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2.1 Matrix Action 



The matrix model that we will study in this paper is given by the gaussian integral over 
a conjugate pair of A/" = 4 supersymmetric matrices Q and Q, acting on a certain finite 
dimensional vector space. To assemble the Af = 4 supermultiplets, we introduce four an- 
ticommuting coordinates ip\ i = 1,..,4 and their hermitian conjugates ipj, and write the 
two matrix variables as superfields Q{ip,ip^) and Q{ip,ip^). We will assume that ip and ip"^ 
satisfy the algebra \^ip\ipj^ = 5*. Besides the anti-commuting coordinates ip\ we now also 
introduce four mutually commuting matrices Z", a = 1, ..,4, of the same size as Q and Q, 
and their hermitian conjugates Z^. We will specify the precise form of these matrices 
momentarily. We combine the matrices and anti-commuting variables as 

= {Z'^\i)') (2.1) 

which can be viewed as a system of coordinates on C^'^. The action for the gaussian matrix 
model now takes the following simple form [12] 

5mm(Q, Q) = Tr {XuQZ'QZ'^ . (2.2) 

Here X/j is a pairing which is anti-symmetric (resp. symmetric) on the bosonic (resp. 
fermionic) part of C^'^, and Tr denotes the trace over the supersymmetric vector space on 
which the matrix superfields act. We will specify this vector space below. 

Now let us specify the four matrices Z". Initially, we introduce the and their 
hermitian conjugates as bosonic oscillators, which satisfy the canonical commutation 
relations \Z"-^ Zj^'^ = 5^. The representation space of this algebra looks like the Hilbert 
space of four simple harmonic oscillators, on which the act like annihilation operators 
and act as creation operators. Summarizing, the oscillator algebra is: 

[Z^Zj]=5^ ; {i^\^}}=S} (2.3) 

which gives an oscillator algebra representation of gl{4:\4, C). This contains ps/(4|4, C), the 
complexified superconformal algebra. We denote by "Hc^k the Fock space of states generated 
by the creation operators Zj. The Hilbert space "^£414 is the linear space spanned by all 
the fuzzy points on a non-commutative C^'"^ [TTIIM]- Each basis state in the Fock space 
represents one Planck cell of the non-commutative space, and since C^'^ is non-compact, 
the associated Hilbert space ^£414 is infinite dimensional. 
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To make the Hilbert space finite dimensional, we will now take the Kahler quotient 

ei4//f/(l) = CF^\\ (2.4) 

where the U{1) acts by uniform phase rotation on all supercoordinates . The projective 
space CP^'^ is compact, and its non-commutative realization has a finite number of Planck 
cells. We should thus expect to find a finite dimensional Hilbert space. The U{1) symmetry, 
that features in the Kahler quotient (12.41) . is generated by the homogeneity operator 

2)o = ^iZ- + ^,V. (2.5) 

This operator has an integer spectrum, given by the sum of the Z-^ oscillator levels. To 
perform the Kahler quotient we consider eigenstates of Hq at some fixed level A^: 

Ho\^) = N\^) (2.6) 

Note that the level constraint Hq = N indeed eliminates one complex dimension: it fixes 
the absolute value of Z"' but also implements the U{1) invariance under phase rotations 
_!. Q^o'Z^ . The condition Hq = N plays the same role as the D-term constraint of the 
usual gauged linear sigma model realization of CP^'^. 

As anticipated, taking the quotient produces a finite dimensional Hilbert space, which 
we will denote by Tii^paii^N). States of TiQ-pa^N) are created by homogeneous degree 
polynomials in the creation operators, acting on the vacuum state. Counting only the four 
bosonic oscillators Z^, this represents a linear space of dimension 

dim-Hcp^iN) = liN + 1){N + 2){N + 3) = kN- (2.7) 
o 

Taking into account the fermionic oscillators, the dimension of "Hj^pSK ( A^) is: 

g 

dimn„3\4{N) = kN + 4fc7v-i + QkN-2 + 4/c7v-3 + kN-4 = -N {N^ + 2) = K^. (2.8) 

We can now specify the form of the matrices Z", by identifying them with the matrix 
elements of the corresponding oscillators between states at some finite level A^. Note, 
however, that Z'^ does not commute with Hq but reduces the level A^ by one. The Z"'s 
thus define maps from 'Hcp3(A^ + 1) to T-Lq^sIN). In other words, the are non-square 
bosonic kj\f x /cat+i matrices. In the supersymmetric case we view the Z's as K]\f x K]\f^i 
matrices. To write the gaussian matrix action (12. 2p . we therefore need to define the matrix 
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variables Q and Q as non-square matrices of size 



Q,Q e Ma.t{kN+iy< kjy). 



(2.9) 



More precisely, this is size of the the lowest superfield component of the matrix variables, 
when acting on the lowest superfield component of the Hilbert space. Since ip'' and xl^j carry 
homogeneity charge —1 and 1, higher superfield components have shifted ranks relative to 
the lowest components. The matrices Q and Q then fill out -R'at+i x Kjsf matrices. In what 
follows we shall often leave the extension to the supersymmetric case implicit, so we write 
all expressions in terms of the bosonic matrix model: 



The variables Q and Q are for the rest arbitrary matrices. A convenient representation of 
the space of arbitrary /cat+i x matrices is as the space of homogeneous polynomials in 
and of degree 1, that is, polynomials in which each term contains one more creation 
operator than annihilation operator, with the relation that = - since any state in 

^cp^(^) is mapped to after acting 1 times with the Z's. Apart from this restriction, 
or after taking the large limit, we can thus view the matrix variables Q and Q as arbitrary 
sections of the degree — 1 hue bundle 0{-l) defined on CP^'^. At finite A^, they are sections 



Finally, let us make a specific choice for the pairing X/j. To this end, we decompose the 
four coordinates into two two-component variables cu" and iTa as 



Tr(4^gz°gz^) 



(2.10) 



of C(-l) defined on fuzzy CP^I^ 





(2.11) 



In this notation, we choose the matrix X/j to be of the following form 




(2.12) 



where rjij is a four index symmetric tensor. It defines a pairing on C' 



'414 



{Z1Z2) = (vriTTa) + [U1U2] + {i^ii^2)- 



(2.13) 
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where we introduced the spinor inner products 



{txitX2) = TXiaTT2b ] [u;iu;2] = f^ab^^Wg. ; = ?7ij?/^}?/;^. (2.14) 

In the following, we will freely raise and lower the spinor indices with the help of the 
corresponding_e symbol. The canonical commutation relations in the two-component spinor 
notation reaqj 

[na,Tx\]=eah ; [uJa,ujl]= e^^. (2.16) 

The eigenvalue condition on the homogeneity operator, that defines the Hilbert space 
T-Lqjps^N), takes the form 

Ho\^) = (vrtTT^ + ooloo' + 4^') \^) = N\^). (2.17) 

Since Hq keeps track of the spinor helicity, it is sometimes also called the helicity operator. 
We see that states in Ti^pa^N) possess a large net helicity equal to A^. 

In the following, we will study the correlation functions of special bi-linear 'current' 
operators computed in the supersymmetric matrix model at level A^. We will focus on the 
leading behavior in the limit of large A^. In this limit, the size of the Planck cells, i.e. the 
scale of non-commutativity, tends to zero relative to the total size of the projective space 
CP^'^. We can thus expect that the large matrix model shares properties with some 
continuum field theory. In the naive continuum limit of the matrix model action fl2.10p . 
the trace over 7iQ-p3{N) turns into an integral over commutative twistor space CP^. The 
resulting free field action takes the form S = J^^s QDQ with bosonic kinetic operator 

o o o 

D = la^Z"^ = Tt"— + COaj^. (2.18) 

dZfS oT" duja 

We should point out, however, that the matrix theory is in fact better defined than the 
continuum field theory with the kinetic operator fl2.18p . As we will see, the continuum 



^There is an unfortunate clash of notation for the square brackets. The square brackets [....] 
around two spinors without a comma in the middle denotes the anti-symmetric pairing of two left 
handed spinors, where the brackets with a comma in the middle denote the usual commutator. 
As another warning to the reader: the notation for hermitian conjugation here contains a raising 
operation for the indices: the hermitian conjugate of tTq and uja is in fact not equal to vrl and w^, 
but rather 

(vrjt = e'^Kl (a;,)t = e-^J^^ (2.15) 
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theory is ultra-local in the sense that the modes only propagate in one direction, and stay 
localized on the other directions. Ultralocal theories do not really exist as continuum the- 
ories, since they typically lead to amplitudes that contain factors proportional to 5(0), the 
Dirac delta-function evaluated at 0. In the matrix theory, this divergence is automatically 
regularized. As we will see, this means that the large N limit of the matrix model has to 
be taken with some care, so as to preserve the UV regulator scale. 

Of course, our motivation for studying the large limit of the matrix model (12. 2p 
is not to regulate some unusual looking ultra-local theory on a complex 3-dimensional 
projective space. Projective superspace CP^'^ is the A/" = 4 supersymmetric version of 
twistor space. The matrix variables Q and Q can thus be viewed as sections of bundles on 
fuzzy twistor space, and via the twistor correspondence, they will then acquire a space-time 
interpretation. 

The twistor correspondence is based on the observation that, given a two component 
spinor tTq and a space-time point x"" on complexified Minkowski space, one can define a 
corresponding two component complex spinor u"" via = ix°'""ira. This relation is invariant 
under simultaneous complex rescaling of the spinors tTq and w". For a given point x, it 
defines a CP^, called the twistor line associated with x. Similarly, a point (x"",6'*") in 
(chiral) Minkowski superspace specifies a bosonic CP^'° in supertwistor space CP^'^, via 

UJ^ = ij' = e'^'TTa. (2.19) 

Based on this space-time correspondence, we may thus expect that a suitable class of 
correlation functions of the large A^ matrix model take on the form of space-time scattering 
amplitudes. In the following, we will consider two situations. In the first case, we prescribe 
that the variables Q and Q, in addition to being kN+i x matrices, also carry an in- 
dex that transforms under the fundamental representation of an internal symmetry group, 
U{Nc). In the strict large A^ limit, the current correlation function then reproduce the 
MHV amplitudes of A/" = 4 SYM theory with gauge group U{Nc). Secondly, we will study 
a specially tuned large A^ scaling limit, where we simultaneously zoom in on a small region 
within the projective superspace, in such a way that the scale of non-commutativity is kept 
fixed. The correlation functions in the resulting double scaled matrix theory reproduce the 
MHV amplitudes of gravity, where the short distance cutoff coincides with the Planck scale. 
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3 Symmetries and Currents 



In this section we will take a first look at D the matrix model kinetic operator which acts 
via: 

DQ = la^Z'^QZ^. (3.1) 

Since each has homogeneity one, and thus changes the level by one, D defines a linear 
map between two spaces of matrices: 

D : Mat(A;7v+i x ^ Mat{kN x kN+i). (3.2) 

We see that the support and image space have the same dimension. Thus we should expect 
the D operator to be invertible as long as the anti-symmetric form is invertible. We 
will investigate the inverse of D later on. 

In this section, we begin with a study of the symmetries of D. As we will see, this 
symmetry group is very large. Via the analogue of Noether's theorem, this implies that 
the matrix model contains a rich collection of current operators. We then study some 
preliminary aspects of correlators built from the symmetry currents of the theory. These 
generate a u{Nc) x gl{k]\;) current algebra. In the later sections, we show that correlators 
of suitably defined currents compute scattering amplitudes. 

3.1 Global Symmetries 

To frame the discussion, we will look at the symmetries of the matrix model through the 
lens of the twistor correspondence f l2.19p . Hence we will view the Z"'s as providing a 
twistor parametrization of space time. The symmetry transformations then acquire the 
interpretation as space-time conformal transformations. 

The main benefit of the twistor parametrization of space-time is that the conformal 
group is generated by linear vector fields. Even in the usual discussions of twistor space, 
it is standard to introduce canonically dual twistor variables Z^, with [Z", Z^] = M^, and 
write the symmetry generators as = ZaZ^. Via the commutators, these manifestly 

generate gl{4:,C), which contains the 4D complexified conformal algebra sZ(4,C). A choice 
of space-time signature amounts to picking an appropriate reality condition. For Euclidean 
signature, one imposes the reality requirement Z^ = Z^- This naturally leads to the 
commutation relation (12. 3 p and the construction of the finite dimensional Hilbert spaces 
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The 15 conformal generators act on 'Hcp3(A^) via the traceless 4x4 matrix of operators 

Mo^p = ZiZp (3.3) 

where we lower indices via la/s so that Za = la/^Z^. We can associate to each conformal 
generator a linear operator that acts on the space of functions $ on the non- commutative 
twistor space (that is, on the space of linear operators $ acting on the Hilbert space 
n„3{N)) via 

Mlp^ = Map<!> - (3.4) 

A simple calculation shows that these operators all commute with the action of D on $: 

[D,MI^]=0. (3.5) 

Based on this equation, it looks as if the kinetic operator D of the matrix model preserves 
the full conformal invariance. However, the matrix model action also involves a trace over 
the Hilbert space 1-Lc.p'-^{N) . In order to be a true symmetry of the action, a charge needs 
to be hermitian with respect to the inner product, that is used in defining the action. 

Hermitian conjugation provides a reality condition, leaving us with the generators m(4) C 
gl{A,C). This is further broken to 5*0(5) by the introduction of the anti-symmetric bitwistor 
Iaf3U The hermitian charges that leave the bitwistor invariant are Ai{a(3) = | , 
which are the 10 generators of S0{5). So these are the true global symmetries of the matrix 
model. As we will see, from the space-time perspective, the matrix model indeed naturally 
lives on the four sphere S*^. 

The space-time interpretation becomes more evident when we write the symmetry gen- 
erators in terms of the two component spinors iTa and u"". We can then distinguish trans- 
lations, conformal boosts, Lorentz rotations and the dilatation generator 

Paa = (^Ua ] Jab = ^l^b) ^ . 

D = -{ulu^ - ny) (3.6) 

Kad = TT^^a ; Jab = '^(a'^b) , 

Hermitian conjugation leaves the SU{2) x SU{2) generators J and J intact, but acts on 
■^The bitwistor transforms as a 6-component vector under so(6) ~ sn(4). 
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the other symmetry generators via 



pL = K"' ; D^ = -D (3.7) 

The hermitian charges that also leave the anti-symmetric pairing (12.121) intact, are the 
5*0(5) generators, J, J and an additional four hermitian generators given by: 

Vaa = Paa + K^" (3.8) 

Thus, only the SO (5) subgroup of the conformal group is unitarily realized. Adding super- 
symmetry is easy. The supersymmetric kinetic operator T>Q = IjjZ^QZ-^ commutes with 
the generators of the complexified super conformal algebra ps/(4|4). The supersymmetric 
hermitian charges that leave V invariant generate the isometrics of the supersymmetric four 
sphere S"''^. 

The symmetry group of the kinetic operator D is in fact much bigger than the global 
isometry group SO (5). Namely, we can consider operators that (similar to D) act on 
operators $ with oscillators from the left and from the right. This allows us to define 
another class of symmetry generators in the form of linear operators Xap that act on 
matrices $ via0 

■ $ = Z[a<f 4] • (3.9) 

One easily verifies that D also commutes with X^/s 

[D,X^^]=0. (3.10) 

The operators Xap, and the fact that they commute with D, will play an important role 
in what follows. For reasons that will become apparent shortly, we will call X"^ position 
operators. Equation (13.101) shows that D is ultra-local, in the sense that it does not shift the 
value of the position operators. This gives us a first precise hint that up to a quantifiable 
amount of Heisenberg uncertainty at small scales, the large N matrix model preserves a 
notion of space-time locality. 

■^There also exists a complex conjugate set of operators X^^^ = ZyJ^Zj^-^, which also commute 
with D. These complex conjugate fields are less relevant for our later discussion. 
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3.2 Currents 



In the accompanying paper [16], we argue that the matrix model appears as part of a 
larger theory that also contains a gauge field A that lives on fuzzy CP'^'^. The fields A 
are linear operators which act on the Hilbert space 'Hj-p3|4(A^). Here we will just focus 
on the subsystem obtained by setting the non-commutative gauge field A to some fixed 
background value, appropriate to the description of scattering states. 

From the perspective of the gaussian matrix model, this coupling is obtained by intro- 
ducing a U{Nc) fiavor symmetry under which the Q transforms in the fundamental and Q 
in the anti-fundamental. In other words each variable Q and Q defines an A^^^- component 
vector of k^^i x matrices. Gauging this symmetry results in the action: 

S = Tr{QD^Q) (3.11) 

where the covariant derivative D_a = D + A acts on fields Q via 

DaQ = 4/3(^° + A")QZf^ (3.12) 

The U{Nc) gauge field A represents some fixed (0, 1) form on the non-commutative twistor 
space. The action (13. lip is invariant under gauge transformations 

6fQ = fQ, 6fQ = -Qf, 6fAc. = [Z^ + /] (3.13) 

where / = /a^Z, Z'^)t^, with an element of the Lie algebra of U (Nc), denotes an infinites- 
imal gauge variation. A basic observation, which will have important consequences later, is 
that even when the color gauge group U{Nc) is abelian, Nc = 1, the gauge transformations 
retain their non-abelian character. Indeed, a U{1) gauge theory on a non-commutative 
space is still non-abelian. As explained in [16], the matrix model enjoys a more general 
gl{k]^) symmetry, which acts by left multiplication on the Q's and right multiplication on 
the Q's. This is automatically gauged due to the presence of the bulk gauge field. So we 
see that there is actually a u{Nc) x gl{kN) gauge symmetry. 

Fixing the background value of A, the coupling between the matrix current and the 
bulk gauge field provides a class of vertex operators for the theory: 

J (A) = Tiil^^A^QZ^Q). (3.14) 

The computation of current correlators then reduces to specifying a choice for the back- 
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ground field A. 

The possible background gauge fields are dictated by the equation of motion ^-"(0,2) = 
J(o,2) which relates the (0, 2) component of the curvature for ^ to a choice of background 
source. Working to linearized order in the bulk equations of motion, we see that there is a 
special class of solutions of the form: 

_ ^ay^ V = ® t"" (3.15) 

where here is an element of the Lie algebra u{Nc) x gl{ki^) and for now, Va is some 
arbitrary function of and Z}^. This choice is an accord with the fact that if we had 
provided a source for the gauge field by activating a vev for Q and Q, the gauge symmetry 
would have been broken. The zero energy configuration would then have been a gauge field 
of the form A"" = Z'^Vr — VlZ". Using the residual generators of gllk^) not contained in 
u{kN), this can be put in the form of equation fl3.15p . 

Plugging the linearized gauge field (I3.15P back into the matrix model action (13. lip , the 
corresponding moment of the current reads 

J{V)=Ti [l^p Z'^Q Z^VQ^ (3.16) 

In this expression, we recognize the kinetic operator D that appears in the free action. So 
we can adopt a more compact notation, and write 

J{V) = Tt{VQDQ). (3.17) 

In the subsequent sections, we will be interested in computing the correlation functions of 
a product of several of these currents in the matrix model. 

3.3 Current Algebra 

Due to the appearance of the kinetic operator in the definition (I3.17p . the currents JT{V) 
vanish on shell, i.e. whenever DQ = 0. As a result, if we study the correlation function 
of a number of current operators, a given current J^{Vi) can be non-vanishing only at the 
location of other current insertions. From the form (I3.17P of the current, we see that J^{V) 
represents the effect of performing "half" a gauge transformation 

5vQ = VQ ; 6vQ = 0. (3.18) 
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So to compute current correlation functions, all we need to do is to perform this substitution 
in the functional integral. 

To make this explicit, consider two neighboring currents. Since the Q and Q are just 
free fields, we can perform a single Wick contraction 

Ti(yiQDQ^ Ti(q{DQ)V^ = Tr (ViQ{DQ)V2^ (3.19) 
I I 

The right-hand side again looks like a current. So we derive the operator product relation 

J{Vi)J{V2) = J{V2Vi) (3.20) 

where the underbracket denotes a single Wick contraction. In analogy with continuum field 
theory, we can thus define a commutator algebra of currents by subtracting the two ways 
of performing the Wick contraction 

[JiV,),JiV2)]=j{[V2,V,]) (3.21) 

The current algebra is isomorphic to the Lie algebra u{Nc) x gl{ki\f) of the V matrices. 

The computation of the correlation functions of currents thus completely trivializes, 
especially when we take the large limit. In this case, we can perform the successive 
single Wick contractions to find 

(^J{Vi) . . . J{Vn)^ = tr(K ...V2Vi)+ permutations (3.22) 

where the sum is over all orderings of the symmetry generators as well as multi-trace 
contributions. So all the physics goes into determining the natural set of V generators that 
we should consider. 

As the reader will have noticed, even when the color gauge group is U{1), the current 
algebra remains non-commutative. Indeed, the U{1) acts not just via phase rotations, 
but also as a diffeomorphism on the fuzzy twistor space. This is a first indication that 
the theory may contain a gravitational sector. To isolate the gravitational physics, it is 
natural to focus on states which in the commutative context would be neutral under the 
U{1). As explained in [TB], this is accomplished by introducing a compensator gauge field 

and viewing the matrix fields Q and Q as bifundamentals under a non-commutative 
u{Nc) X n(l) gauge symmetry. The covariant derivative of equation f l3.12p is then replaced 
by D^^Q = Iai3{Z°'+A°')Q{Z^+A^). The gravitational gauge symmetry then corresponds 
to the linear combination of m(1)'s which acts on the defects via the adjoint action. Note 
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that in the commutative context, the mode would have been neutral under this adjoint u{l) 
action. In our context, this acts via pure diffeomorphisms. 

Let us next consider the corresponding gravitational currents. There are two glik^) 
symmetry currents, given by A'^ = Z'^V and = VZ^, where V and V are functions 
of Z and The adjoint action corresponds to setting V = —V. This leads to an 
additional set of gravitational currents: 

r{V)= Tt[[V,Q]DQ^ (3.23) 

The analogue of equation f l3.24p is then: 

T{Vi)T {V2) = T{[V2,Vi]) (3.24) 

In this case, the V^s directly act via commutators. The commutator algebra of the V^s 
should be viewed in the commutative limit as the algebra of vector fields on twistor space. 
We will make this more precise when we turn to the computation of scattering amplitudes. 
To this end, we now turn to the space-time interpretation of the matrix model. 

4 Space-Time 

Having shown that the gaussian matrix model enjoys a number of symmetries, in this section 
we turn to their 4D space-time interpretation. In particular, we determine a fuzzy twistor 
correspondence between points of a 4D space-time and fuzzy CP^'s. Using this interpreta- 
tion, we can view the matrix variables Q and Q as fields on non-commutative twistor space. 
Via the twistor correspondence, they will then acquire a space-time interpretation. 

4.1 Twistor Lines 

Given the appearance of the symmetry algebra 5*0(5), we should expect some connection 
with space-time physics. Here we develop the notion of a "coherent state" |x. A) which is 
associated with a spacetime point x"" and a local coordinate A on a CP^. The extension to 
the supersymmetric situation will be straightforward, and is discussed in p^fTT] . To begin, 
we start with a normalized state |0,0) which up to an overall normalization is the unique 
state annihilated by the oscillators cj" and 712- Acting by S0{5) generators, we can sweep 
out the rest of Hpj. The states of a fuzzy CP"^ are obtained by acting with J, the su{2) 
subalgebra built from just the vr oscillators. We refer to a holomorphic point on this CP^ 
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as a state |0, A) which satisfies: 

e^'^XaTib |0, A) = and co^ |0, A) = . (4.1) 

where (Ai,A2) are homogeneous coordinates of the commutative CP^ and A = A2/A1 is 
an affine coordinate. The space of all |0, A)'s are mapped to each other via the 5*0(4) 
generators J and J. The equivalence class of all such states is then a fixed point of 5*0(4), 
corresponding to the south pole of an 5*^. 

Starting from the south pole of the 5^, we can now sweep out the remaining states by 
5*0(5) generators. Acting via x ■ V = x"'°-Paa + XaaK"-"' of equation (13.81) . we obtain states: 

\x,X) =exp{ix-V)\0,\) (4.2) 

In this way we build up a spin N/2 su{2) bundle fibered over 5^. Note that the transfor- 
mations exp{ix ■ V) are unitary, and do not alter the norms of states. 

The fiat space limit corresponds to the Wigner-Inonii contraction of 50(5) where we 
rescale the generator Paa relative to K"'"-. In this limit V ^ P, and the states \x, A) satisfy: 

e'^%7r,|x,A) = (4.3) 
{00^ -ix^''7ra)\x,X) = (4.4) 

The second line is nothing but the usual twistor equation associated with a space-time point 
but now interpreted as a holomorphic operator equation. In other words, we recover 
the expected correspondence between a point x°'°' of complexified Minkowski space and a 
(fuzzy) CP^ 

By a similar token we can introduce bra states (x, A| = (0, A| exp{—ix ■ V). We provide 
the precise definition of (0, A| in section [51 In the fiat space limit we obtain bra states 
annihilated by ul — iXaaT^"^"". Note that both the bra and ket states can be extended to 
holomorphic x"""". This is an important feature of twistor geometry which is preserved by 
the matrix geometry. 

Having established a connection with classical twistors and the 4D continuum space- 
time, let us now discuss some additional features of commutative twistors. See for 
additional review. In twistor theory [211I2S], the identification between space-time points 
and complex lines in twistor space is a correspondence at the level of holomorphic geometry. 
Complexified conformally compactified Minkowski space is given by the zero locus of the 
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Klein quadric in CP^ 

e^p.sX'^^X^' = . (4.5) 

Here X""^ = — is a four index anti-symmetric tensor, defining the six homogeneous 
coordinates of CP^. The constraint (14. 5 p is automatically solved by introducing a pair of 
points in twistor space, with homogeneous coordinates U"' and V^, via 

Since two twistor points U and V determine a line Z = all + bV in CP^, one recovers the 
map between space time points and twistor lines. 

The homogeneous coordinates X"^ are sensitive only to the conformal structure of space- 
time. Conformal symmetry is broken by designating a choice of two index anti-symmetric 
bitwistor, called the infinity twistor, denoted by Jq^ = —Ipa- The 'inverse' bitwistor is 
denoted by I^^ = ^e°'^'^^Iys- The infinity twistor defines an anti-symmetric pairing 

{ZW) = lapZ'^W^ (4.7) 

and allows us to raise and lower the index of the twistor coordinates Z" via Z^ = 
With the help of the infinity twistor, we can define affine space time coordinates 

Xe^p = Xo = r^X^^ (4.8) 

The matrix Ijj that features in the gaussian matrix model action is the infinity twistor of 
S*^'^, the supersymmetric four-sphere. 

We have already encountered this formulation of the twistor correspondence in the non- 
commutative setting, in the form of the operator^ X^/^-^ = Z^a^Z'^^-^ introduced in equation 
(13. 9p . We now see the space-time significance of these operators: they allow us to define 
the notion of operators $(x) that are localized at a given space-time point x. In analogy 
with (14. 8p . we shall sometimes refer to an eigenstate of X as a matrix $(x) which satisfies: 

X"^.<|.(x) = x"^X°-<l>(x) (4.9) 

where Xq = I^-^X^p and where Xap are c-numbers. Alternatively, we could have defined 
local operators ^(x) as operators that satisfy the space-time coherent state conditions from 

^Here we temporarily put a " on the coordinate operators X^p, to distinguish them from 
ordinary c- number space-time coordinates. 
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the left, and the hermitian conjugate conditions from the right 

{u^-ix''^7ra)^ix) =0 ; (^{x){u^^ - ix^'^Trl) = (4.10) 

It is not difficult to show that the two definitions (14. 9 p and (I4.10p of operators $(x), that 
are localized in space time, are equivalent]^ The precise relation between the eigenvalues 
Xai3 in (14. 9 p and the fiat space-time coordinates x°'°' that appear in the twistor line equation 
(l4.1Up is given in Appendix A. 

The commutator of Xaa with the 5*0(5) rotation generators Vaa = Pda + K"''^ reads 

[^da,n6] =ea6e,^Xo. (4.11) 

Near the south pole region, where Xq is maximal, we can approximate Xq by its maximal 
eigenvalue. After rescaling Xaa to Xaa = Xaa/Xo, this relation yields the Heisenberg com- 
mutation relation between momenta and coordinates. Finally, since the position operators 
Xai3 commute with the kinetic operator D of the matrix model 

[D,X^p] = 0. (4.12) 

Hence, the i5 operator maps the local operators $(x) (defined via the eigenvalue equation 
(14. 9 p or equivalently, the coherent state condition (I4.10p ) to another local operator (D$)(x) 
defined at the same space-time point x. The kinetic operator D thus acts along the twistor 
lines. 



4.2 Planck Scale 

We have seen that the non-commutative theory allows for the introduction of position 
operators with a continuous spectrum of eigenstates. Of course, this does not mean 
that the matrix model defines an exact local theory. Indeed, from the perspective of the 
4D space-time, the fuzzy twistor space corresponds to truncating the angular momentum 
on the S^. This limits the angular resolution of the 4D theory. The number of independent 
spherical harmonics on 5^ at level N is of order N^/12, with corresponding resolution area 

^The proof is the same for the left twistor Hne equation and its conjugate, and goes as follows: 
= e^p^sZ^X^^Hx) = {e^p^sZ^x''^)X^<^{x) 
Since X^ commutes with Z^, and is invertible, this imphes that [eaiB^ysZl^ x^^)^{x) = 
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Another way to see the presence of this minimal length scale is by evaluating the overlap 
of position eigenstates (f/| = (x, A| and \ V) = \y,^), for different space-time points x and y: 



= (x,A|y,0 = (0,A|e-^^■V^■^|0,O• 



(4.13) 



The composition q-^^-'T' ^w-v jg again an 5*0(5) rotation. At small displacements, it corre- 
sponds to a translation in the direction r^a = {x — y)da- Next, we can use the fact that any 
S0{5) rotation operator can be factorized as a product -R5 = {R4)Rq{R/^) where each {R4) 
factor is an 5*0(4) rotation, and where Rg is the special rotation matrix Rg = e^^^'^ with 
Vaa a unit 2x2 matrix proportional to Vaa- The diagonal 5*0(5) rotation Rg can be thought 
of as the rotation that transports the point x along a great circle to y. The rotation angle 
is the arc length 

e=\x-y\/i (4.14) 
Rg acts on the four twistor coordinates via the simple rotation (see for example [37] ) 



Rg 













— > 







/ 



/ COS I TTi+i sin I c^i \ 
cos I 7i2—i sin | UJ2 
cos I Ui+i sin | tti 
\ cos I — i sin 1 / 



(4.15) 



Ignoring for now the 5*0(4) part of the rotation, one can easily compute the matrix element 
by letting this transformation act on the ket state |0, A). Using that both the bra and ket 
state contain only it oscillators, one immediately finds that the answer collapses to 



(0,A|i?,|0,O 



cos 2 



N 



(0,A|0,O 



(4.16) 



This is the expected behavior of an 5*0(3) transformation with rotation angle 6 acting on 
a spin A^/2 representation. 

We are interested in the leading behavior at large and small 6: 

N\x - y\'^- 



cos 



\x~y\ 
2i 



N 



— )■ exp 



(4.17) 



Here we introduced the UV length scale ipi via 



N 



(4.18) 



The parameter ipi represents a short distance cutoff for our theory. In the last step in fl4.17p 
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we took the large A^, large ^ limit while keeping Ipi very small but finite. Finally, we can 
verify that the 5*0(4) part of the 5*0(5) rotation indeed drops out, because (i) the states 
|0, .^) and (0, A| do not depend on the u oscillators, and (ii) the delta function enforces that 
the rotation parameter vanishes anyhow. 



5 Chiral Field on CP^ 



Having presented the gaussian matrix model and established that it retains a natural 4D 
space-time interpretation, we would like to study the correspondence between correlators 
of the matrix model and 4D physics. To this end, in this section we compute the exact 
form of the propagator for a chiral boson or fermion system on a fuzzy CP^. 

On a commutative CP^, the chiral free field action takes the form 

S= I 4>d(j) (5.1) 

The chiral fields (p and cf) can either both be fermions or bosons. In principle, they can carry 
arbitrary half integer spin s and 1 — s, respectively. In the following we will mostly restrict 
to the spin half case s = |, so that Q and Q are both sections of the degree —1 bundle 
0{—l). Using projective coordinates (7ra,vfb), a = 1,2, on CP"*^, the d operator reads 

d = vr A (5.2) 

diTa 

In this section, we are interested in constructing the non-commutative version of the prop- 
agator of the chiral fields. In other words, we will be looking for the analogue of the 
Green's function A(7r,A) associated with the d operator. Using the projective notation 
(ttA) = e'^^AaVTb, for the difference between two points, our task is to solve the equation 

M(7r,A) = 5((7rA)) (5.3) 

In the commutative theory, this is trivially solved vi 

A{n,X) = j\r- (5.4) 
(ttA) 

As we will see, although finding the non-commutative analogue of this expression takes a 
bit more work, the end result will be almost as simple. 



'^In affine coordmates = 7^2/'^! and A = A2/A1, it takes the more familiar form A(^, A) 
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Non-commutative CP"^ is described by oscillators (7ra,7rJ), with a = 1,2, satisfying the 
canonical commutation relation [na, vr^] = eab- The oscillators act on finite N+1 dimensional 
Hilbert spaces 'H£pi(A^), specified by the level constraint 

e'^V>b|^)=iV|^). (5.5) 

As before, we can think of the Hilbert space 'H^pi{N) as the space of points on the CP^ 
The chiral fields Q and Q represent arbitrary homogenous polynomials in the creation and 
annihilation operators tTq and vrj of a given degree specified by their spin. In the spin 1/2 
case, they are taken to be homogeneous functions with one more vr''' than vr. Hence, the 
fields do not act within the same finite Hilbert space: and both act as linear maps from 
'Hcpi(A^) to Hqpi^N+I), and can thus be viewed as arbitrary (A^ + 2) x (A^ + 1) matrices. 

The non- commutative version of the action (15.11) reads 

5 = Tr(0e"V,</.7rb), (5.6) 

where the trace is taken over 'Hcpi(A^+l). Note that the kinetic operator 

d = e^^'n^Lnn (5.7) 

via the action of the commutator on the tt^ dependence of 0, indeed defines a direct analogue 
of the Dolbeault operator (15. 2p . See [381 - HO] for further discussion on the form of the fuzzy 
Dolbeault operator. 

Given the action, we can start to compute correlation functions. For this we need 
the explicit form of the propagator A = 1/d. Mathematically, the 9-operator defines a 
linear map from the space of (A^ + 2) x (A^ + 1) matrices to the space of (A^ + 1) x (A^ + 2) 
matrices. Since the support and image have the same dimension, this map is expected to be 
invertible. This is indeed obvious from the oscillator representation: d has no zero modes, 
since always contains at least one tt"^ oscillator. We can thus define the non-commutative 
version of the propagator as the inverse of this linear map. 

5.1 AfRne Coordinates 

While the commutative theory (15. ip enjoys full conformal invariance, the non-commutative 
deformation breaks the conformal group to the group of global SU (2) rotations acting on 
the doublet of oscillators na- The non-commutative CP^ is indeed equivalent to a fuzzy two- 
sphere. In the following, however, we will not use this global SU{2) perspective, because 
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we wish to preserve the holomorphic properties of the theory as much as possible. To this 
end, we will choose to work in a local affine patch with coordinate A = A2/A1. As we will 
see, this will allow a formulation in which conformal symmetry will naturally re-emerge 
once we take the large limit. 

Let us introduce the following number basis of 'H^pi{N) and its dual 

I") - ' ('^_„i,,;! |0>. (n| = {0|U,)>,)"-". (5.8) 

The normalization factors are convenient for our present discussion. This basis is canonically 
normalized 

(n|m) = 6n,m (5.9) 

where n and m both run from to A^. The Hilbert space 'H^pi{N) contains a continuous 
family of coherent states, labeled by points A on the commutative CP^, defined via 

N 

|A) = ^,(A) 5^A"|n), (5.10) 

n=0 

where 0^{X) is a normalization factor. A geometrically natural requirement is that A is 
invariant under simultaneous transformation tti ^ tt2 and A O A~^. This leads to 

(^nW= ^_In+i (5-11) 
In the large A^ limit, this becomes a step function: 

1 for |A| < 1 
^(A) = { (5.12) 
for |A| > 1 

The states |A) satisfy the coherent state condition 

(7r2-A7ri)|A) = (5.13) 

which shows that A can be thought of as the classical value of the affine coordinate A2/A1. 
We may write (15.131) in a slightly more covariant notation as 

e'^%7r,|A) = (5.14) 

with Aa = (1, A). We can call the states | A) 'position eigenstates', although there obviously 



23 



does not exist any unitary position operator of which they are eigenstates. The state |0) 
corresponds to the position state at the origin A = 0, while the state |A^) corresponds to 
the point at infinity 

|A^) = |oo). (5.15) 

We will sometimes call |0) the south pole state, and |oo) the north pole state. 

At this point it is useful to introduce the non-commutative notion of the affine coordinate 
chart. The main advantage is that the chiral fields will become square matrices. In our 
setting, specifying an affine coordinate system amounts to picking a 'canonical' embedding 
of 'Hcpi(A^) inside of TiQ-pi^N + 1), or equivalently, a projection from Tii^pi^N + 1) onto 
'Hcpi(A^). Choosing the coordinate A = 1^2/ t^i amounts to identifying the states in both 
spaces via the action of the vri oscillator. In particular, position eigenstates are related via 

= vri|A),^,. (5.16) 

This map projects out the north pole state (15.151) . since 7ri|oo) = 0. We will call this coor- 
dinate chart the south pole patch. The restricted Hilbert space, with the north pole state 
projected out, will be denoted by l-i'^^^{N+l). The map (I5.16P provides an isomorphism 

n„.{N)::^n'„.{N+l) (5.17) 

The chiral free fields in the affine coordinate patch are defined as $ = 7ri0 and $ = 7ri0. 
The redefined fields both act as linear maps from T-Li^pi^N) to itself, and thus specify square 
(A^ + 1) X (A^ + 1) matrices. We will use the isomorphism (15.161) repeatedly in what follows. 

As one would expect, the dual Hilbert space is naturally viewed as describing the op- 
posite patch with affine coordinate ^ = tti/tt^. We will call this the north pole patch. In 
adhering to the usual notions of twistor theory, we seek a suitable holomorphic notion of a 
bra state. At first sight, however (since bra states cannot be annihilated by a linear com- 
bination of annihilation operators) there is no obvious dual basis of position eigenstates, 
which are annihilated by the holomorphic operator tti — ^712. We can still define coherent 
states ((^1 via 

N 
n=0 

with 6^{C,) defined in (15. lip . A straightforward calculation shows that the holomorphic 
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coherent state condition is violated at the north and south pole 



(eKvTi -e^2) = ^.(0(oi - ^.(r')(oo| (5.19) 

Here (0| and (oo| = (A^| denote the dual north and south pole state in 'H*pi(A^ + 1). The 
dual north pole state (0| is located at ^ = and the south pole state (A^| = (oo| is the place 
where .^ = 00. Both states will play a special role in what follows. Note that the inner 
product pairs the dual north and south pole states with their polar opposites 

(0|0) = (cx)|cx)) = 1, (oo|0) = (0|oo) = 0. (5.20) 

In this sense, our inner product is similar to the BPZ inner product in the radial quantized 
formulation of 2D conformal field theory. 

In the large N limit, the factor 6*^(0 becomes a step function: it is equal to 1 on the 
northern hemisphere where |^| < 1, and vanishes on the southern hemisphere where |^| > 1. 
So after taking the large hmit, eqn fl5.19p reduces to 

f (0| for|e|<l 
m7ri-^7r,)={ (5.21) 
[ (cx)| for 1^1 > 1 

This is our desired intermediate result. It shows that the dual coherent state (.^| are position 
eigenstates, modulo a source term localized at the corresponding pole. 

Let us compute the overlap between the position eigenstates. A direct calculation shows 
that 

This equation reveals, as expected, that ^ and A are reciprocal affine coordinates. Note, 
however, that the pole in the denominator is spurious. The step functions do not allow ^ 
and A to be located on the same hemisphere: whenever they do, the numerator vanishes. 

The standard way to overcome this obstacle is via analytic continuation. Consider the 
overlap (^i|A2) and let Ai be the reciprocal coordinate to .^1 = Aj~^. We wish to define the 
south patch state (Ai| via analytic continuation of the north patch state (^i|. However, 
here we meet a subtlety. At infinite A^, the step functions are non-analytic at the equator, 

^Here we use the action of -K2 to define a canonical embedding of the dual Hilbert space ^^pi (A^) 
inside 1W* (A^+l). 
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while at finite the step functions On{C,) are perfectly analytic. Our approach is to first 
take the large TV limit, and then analytically continue. See Appendix B for more discussion 
of this issue. 

The upshot is this: we define the state (Ai| such that its overlap with IA2) is given by 
the analytic continuation of f l5.22p . starting from the region 6{^) — 6{X^^) = 1. Hence in 
the strict large N limit, we hav^j 

Moreover, via the analytic continuation of fl5.2ip . we learn that (A|, as defined this way, 
solves the ket state condition (A|e'^^Aa7r;, = (0| up to terms which vanish in the large 
limit. This completes our construction of the state (A|. 



5.2 CP^ Propagator 

The construction of the propagator of the chiral free fields is now almost as simple as in 
the commutative case, or possibly even simpler. We first need to define the notion of a 
holomorphic delta function. Let A be a point on the commutative CP^, with associated 
coherent state |A). Our proposed definition for the projective delta function localized at 
the point A is as follows 

<5((7rA)) = |A)(0|. (5.24) 

Let us motivate this definition. Eqn f l5.24p defines a projection onto the position eigenstate 
at A, which is similar to how a commutative holomorphic delta function acts on the space of 
functions. The right-hand side explicitly involves the special state (0| which represents the 
point at infinity of the affine chart tti = L The projective delta-function on the left-hand 
side seemingly does not depend on such a choice - but of course it does once we choose an 
affine chart on CP^. 

Given this definition of the projective delta function and the result fl5.19p . we now have 
a natural candidate for the propagator 

A(7r,A) = |A)(A| (5.25) 

^In going from the patch near the north pole to the patch near the south pole, we have used 
the fact that the bra states transform as 1/2-difFerentials, and so in passing from one patch to the 
other, transform as (^| — t- ^(A|. 
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The verification is trivial: 

e'^VjA)(A|7r, = |A)(A|e'^''A,7r, = |A)(0| (5.26) 

up to corrections which are exponentially suppressed at large A^. Given the proposed 
identifications, this calculation provides the non-commutative version of eqn (15.31) that 
defines the Green's function of the d operator. 

6 Scattering Amplitudes 

Having studied correlators of the CP^ system, we now turn to correlators of the full gaussian 
matrix model. In this section we propose a direct correspondence between correlators of the 
matrix model and amplitudes of the 4D space-time theory. This correspondence is defined 
in a double scaling limit where we zoom in on a small neighborhood near the south pole of 
the S^: 

N ^oo, £ ^ oo, 4; = ]^ fi^ed (6.1) 

Scattering in the 4D theory proceeds as follows. We prepare states "at infinity" , correspond- 
ing to the boundary of the small patch near the south pole. The rescaled patch defines 
our 4D spacetime for the scattering experiment. To have a notion of lightlike momenta, we 
compute the values of the correlators in Euclidean signature, and then analytically continue 
to lightlike values of the complexified momenta 

Pad = AaAd (6.2) 

for complex spinors A^ and A^. We can then speak of a matrix model current J7i for a 
massless state with a specified momentum pi. The basic dictionary is that a scattering 
amplitude is represented as a correlator of currents in the matrix model: 

i-Ml,...,n = (Jl-Jn) ■ (6.3) 
\ /MM 

Here the correlator is evaluated by performing the matrix integral while taking the double 
scaling limit (16. ip . 

This section is organized as follows. First, we begin with a discussion of the fiat space 
limit, and in particular, how to pass from an abstract correlator of the matrix model to a 
scattering amplitude. To evaluate such correlation functions, we need to construct the CP'^ 
propagator. As we will see, our detailed study of the CP^ example will give a good return 
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Figure 1: Starting from a round S*^, the flat space limit is obtained by zooming in on a small 
patch near the south pole of the geometry, depicted by the shaded red region. This is then 
rescaled, yielding M^. As depicted in the right panel, scattering amplitudes are computed by 
analytically continuing the correlator in the fiat space limit to general complex momenta. 

of investment. Then we construct the asymptotic wave functions. Finally, as a warmup for 
our discussion of graviton amplitudes, we discuss how the model reproduces MHV gluon 
amplitudes. 

6.1 Flat Space Limit 

In order to compute scattering amplitudes, we need to pass to a 4D theory on fiat space- 
time via the double scaling limit (16. ip . In this subsection we discuss in more detail how to 
treat this limit. See figure [1] for a depiction. 

When we zoom in on the region near the south pole region where becomes small, 
or equivalently, where u°' is much smaller than tt^, the S*^ curvature becomes negligible 
and the space-time enjoys an effective translation invariance. The hermitian translation 
operators are given by the generators V of 5*0(5). In the fiat space limit, they are related 
to the generators P and K via the Wigner-Inonii contraction of the so(5) algebra: 

V^P + r'^K (6.4) 

In this limit, V remains hermitian, provided the dagger of P is now P^ = i~'^K. In terms 
of the oscillators, we rescale the u oscillators relative to the it oscillators, while keeping the 
5*0(4) subalgebra invariant. This is also reflected in the equation u"" = ix°'""iTa, so that when 
X has dimensions of length, u is scaled relative to vr. Since the conformal boost generator 
leaves the south pole of the 5*^ fixed, its effect becomes negligible in the scaling limit ( 16. ip . 
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To make the connection with ordinary translations more exphcit, consider the commu- 
tator between the Vaa generator with the space-time coordinate operators Xaa introduced 
in Section HI 

[Xaa.Vl^^] = eabe^^Xo. (6.5) 

The function space near the south pole region is given by linear combinations of eigenstates 
of Xaa/Xo, as defined in ( I4.9p . with eigenvalue Xaa ^1- In this region, Xq attains its 
maximal value, and can be treated like a c- number constant. Hence in the scaling limit, 
"Pj,^ acts like a translation operator. 

Evaluating correlation functions in this limit involves the insertion of the projection 
operator 1m into the definition of the amplitudes: 

1m = / d'^x^\x,p){x,p\ (6.6) 
p 

where the domain of integration for x is over the Minkowski patch near the south pole. 
Here, the states \x, p) are obtained by starting from the CP^ at the origin, and sweeping 
out by the so(5) generator exp(ix ■ V). Similar considerations hold for (x,p|. For the 
most part, such insertions can be ignored. However, when we turn to a discussion of MHV 
graviton scattering in section [7] where the initial states themselves disturb the location of 
the patch (as they are infinitesimal diffeomorphisms) , additional care must be taken. 

6.2 CP^ Propagator 

The D kinetic operator (13. ip of the twistor matrix model essentially reduces to the CP^ 
kinetic operator acting along the twistor lines. This fact can be anticipated by taking the 
naive commutative limit of the matrix model kinetic operator (I2.18p . This operator acts as 
a one-dimensional d derivative along twistor lines. Hence we expect the continuum limit of 
the propagator to be delta- function localized along the directions transverse to this line. 

The propagator satisfies the inhomogeneous wave equation with a delta function source. 
Following Penrose, we pick this delta function source via the pull back to the correspondence 
space. Let 

5{Z, U) = 6\uj^ - 2x"X) SiiirX)) (6.7) 
be the delta function that localizes Z = (w", tt^) at U = (x^^A^, A^). The continuum CP^ 
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propagator is defined as the Green's function tliat solves 



DA{Z,U) = S%Z;U) (6.8) 

It is easily verified that the solution reduces to a projection operator onto the twistor line, 
times the CP^ propagator on this line. 

A(Z.^)^ ^^'-;-f-°) (6.9) 

(ttA) 

We now translate this to the non-commutative setting. 

The notion of the correspondence space relies on the use of complexified space-time 
coordinates, where translations are generated by the operators P rather than their hermitian 
counterparts V. In the following, however, we will be interested in the limit in which the 

gets very large, that is, we zoom in on a small region near the south pole, which in the 
large radius limit approaches fiat space. In this region, the violation of 5*0(5) symmetry is 
minimal. Conversely, the 5*0(5) generators act to a very good approximation as translation 
generators of the Poincare group, which are compatible with the holomorphic data and do 
preserve the form of the Green's function (16. 9p . 

Our construction of the Green's function is modeled after the one employed for the CP^ 
case. We will use the correspondence space parametrization and choose an affine coordinate 
patch X = -^2/111. To every point labeled by (x. A), we can associate a coherent state via 
(14.31) and (14. 4p . Our strategy is to first find the delta function and Green's function that 
are localized at the twistor line at the origin x = 0. We will then find the general solution 
by acting with the 5*0(5) symmetry generators. 

By analogy with the CP^ case, we identify the holomorphic delta function that localizes 
on a point A on the twistor line at x = as S{Z; 0, A) = |0, A)(0, 0|. Here (0, 0| is the dual 
north pole state on the CP^ at the origiiF"! 

(o,o| = £;,MoK (6.11) 

Note that both the state |0,A) and (0,0| are just made up from the vr oscillators. So we 



Here we included a factor of so that, in the large N limit, the overlap of (0,0| with 
south pole states |0, x) at other locations is properly delta- function normalized: 

(0,0|x,0) = (5'^(x). (6.10) 
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can view both states as part of the CP^ Hilbert space at the origin. We can thus carry over 
the results of the previous section, and derive that the Green's function, that satisfies (16. 8p 
with the delta function 6{Z;0,X) = |0,A)(0,0|, is given by A(Z; 0, A) = |0,A)(0,A|. Here 
(0, A I is the state that satisfies: 

(0,A|e'^%7rfe = (0,0| (6.12) 

up to terms which are exponentially small at large A^. To move away from the origin, we 
act by X -V: 

|a;,A) = e^^-^|0,A) ; (x, A| = (0, Ale""-^. (6.13) 

where we assume that the magnitude of x is very small in the flat space limit. 

Following the by now familiar pattern, our definition of the delta function (16. 7p is 



S^{Z;x,X) = \x,X){x,0\ (6.14) 

The right hand side projects onto the position eigenstate \x, A), as the delta function should; 
the bra state (x, 0| corresponds to the north pole of the S"^ associated with the affine 
coordinate system on the twistor line for x. In the fiat space limit, the state (x, A| satisfies 
the inhomogeneous holomorphic coherent state condition: 

(x, A|e"^A,(7rfc - iXabCo'') = (x, 0| (6.15) 

up to terms which are exponentially suppressed at large A^. How unique is a solution to 
this equation? We notice that this equation only involves one linear combination of the 
oscillators. Since there are several oscillators, it might look like this single condition does 
not uniquely fix the state. However, suppose we had found another state that solves eqn 
(I6.15p . Taking the difference with our solution for (x, A| would yield a bra state that is 
annihilated by a linear combination of annihilation operators. Clearly no such state exists. 
Hence our solution is unique up to small correction terms. 

The propagator that satisfies (16. Sp . with the above identification of the delta function, 
is now immediately found to be 

A(Z;x,A) = |x,A)(x,A| . (6.16) 
Verification of the Green's function property follows immediately from the fact that [D, V] = 



31 



0. Indeed, since ^(10, A)(0, A|) = |0,A)(0,0|, we obtain: 

D{\x, A)(x, A|) = D(e^^-^|0, A)(0, Ale"'^'^) = |x, A)(x, 0| (6.17) 

The final expression (16.161) for the CP'^ propagator will be used repeatedly in the following 
sections for the computation of scattering amplitudes. Because D is an invertible map on 
this basis of matrices, we can also invert both this map and the action by exp{ix ■ V). This 
establishes the uniqueness of the Green's function solution. 

A last piece of information we need is the inner product between the special bra states 
{x,\\ with a position eigenstates. Using the earlier calculation of the overlap of position 
eigenstates, we find that 

(xi,Ai|x2,A2) = 5\xu) (A1IA2) = (6.18) 

which should be compared with the continuum version (EH) of the CP^ Green's function. 

The generalization to supertwistor space is straightforward. Starting from the CP^ at 
X = ^ = 0, we have the delta function on the CP^, |0, 0, A)(0, 0, 0|. At small x and 9, this 
corresponds to a point on the supercorrespondence space U = (x'^'^Aq, 0^°'Xa, A^). In the flat 
space limit, the state (x, ^, A| satisfies: 

(x, 6, Xle'^'Xain - ixabu^ - r]ijtlj%) = (x, 6, 0| (6.19) 

Acting by a symmetry generator of 5"^'^, we can move out to a general value of x and 9. 
This symmetry generator commutes with the supersymmetric kinetic operator V, so again 
the verification of the Green's function property is trivial. The propagator is then given by 

{x^,9,,X,\x2,e2,X2) = ^-^^^ (6.20) 
where the delta function is over the A/" = 4 superspace. 



6.3 Space-Time Currents 

The next step in assembling the ingredients of the S-matrix is to determine a physically nat- 
ural basis of currents. These are specified by a choice of background gauge field A"" = Z°'V 
for some V = 'Va® and a generator of u{Nc) x glik^)- The generators u{Nc) define 
currents on color space, while the gl{k]^) generators are deformations of the geometry itself. 



32 



At a heuristic level, we are interested in taking V = V{Z) to be a "locally holomorphic" 
function of just the Z's in the sense that they commute with the holomorphic coordinates. 
In this sense, such functions do not disturb the holomorphic geometry of twistor space. 

Strictly speaking, this cannot really be done on a finite size S*^, and in particular in the 
finite theory. The reason is that all matrices we write down will be a power series in both 
Z and Z'^ , so all currents will inevitably distort the geometry. Indeed, it is precisely this 
feature which suggests a connection with gravity. This is closely related to the presentation 
of the position eigenstates \x,X) in the flat space limit. Recall that in this limit, |x. A) 
is obtained by starting from the south pole state |0,A) and applying exp{ix ■ V). This 
intrinsically links this collection of states to a small neighborhood in the vicinity of the 
south pole. On the finite size 5*^, we could have alternatively started from the north pole 
and rotated by a different 5*0(5) rotation to reach the same point on the S"^. Note, however, 
this operation would not have been holomorphic in the original x"'^, as it involves transport 
from the point at infinity. Hence, when we work at finite A^, the most we can hope for 
is an approximate notion of holomorphy in the V{Z) which becomes exact in the large A^ 
limit. When we turn to a discussion of plane wave solutions, we shall give a more precise 
characterization of such "locally holomorphic" V^'s. 

In the following, we will distinguish two special classes of V generators. The first class 
are the closest analogue of local color gauge rotations acting on Q. 

V^{Z)=T^V{Z) (6.21) 

We can call these transformations local color rotations, because they act on the color index 
of Q but otherwise commute with the holomorphic coordinates (at least locally). They 
therefore do not induce any coordinate shift of the holomorphic coordinates. As we will 
see, the correlation functions of currents associated with this class of transformations will 
correspond to gauge theory MHV amplitudes. 

A second special class of transformations are those that leave the color index unchanged, 
but act non-trivially on the holomorphic coordinates by means of an infinitesimal gl{kN) 
transformation. A natural class of generators are |16j : 

VUZ) = VaaV{Z) (6.22) 

where Vaa = Pda + K"'"' is an so(5) generator. The lefthand side Vaa{Z) is a glikjq) generator 
which contains a single oscillator. Via the commutator, it describes a holomorphic vector 
field on the non-commutative CP'^. Correlation functions of currents associated with this 
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class of transformations will correspond to gravity MHV amplitudes. 

In the supersymmetric case there are additional transformations and associated cur- 
rents. These are given by the purely fermionic sm(4) generators, as well as mixed bosonic 
and fermionic currents. The former can be identified with the gauged R-symmetry of a su- 
pergravity theory with J\f = 4 supersymmetry, while the fermionic components correspond 
to the gravitinos. 

6.3.1 Plane Waves 

To complete our characterization of the y's, we now construct operators corresponding to 
asymptotic states with specified complexified momentum pad = ^a^d, as appropriate for a 
discussion of scattering theory in twistor space. 

The construction of the solutions in the fiat space limit is obtained by viewing \x, X){x, A| 
as a designated projection to a point of the correspondence space. This is of course in accord 
with the identification of the flat space limit projection matrix 1m of equation (16. 6p . A 
momentum eigenstate is then given by summing over the continuum position x of the 
operator e*^'^|x, A)(x, A| in this small patch. 

We now formalize the algebraic conditions for a matrix to be a momentum eigenstate. 
The first requirement is that the propagating mode V is specified by an asymptotic source 
/ via lapA'^Z'^ = f, for f a. kjq x kjq^2 niatrix. Dropping all group theory indices, the 
identification A"" = Z'^V is summarized by the condition: 

DV = f. (6.23) 

Here, / is treated as an a priori arbitrary source for the wave function V. 

To construct a momentum eigenstate, we now further restrict attention to operators 
V{p) with a specified momentum pda- In the fiat space limit, this is designated by the 
condition: 

[Pda,Vip)]=PdaV{p) (6.24) 

where VOp) is viewed as a state in the large limit of the adjoint representation of 
gl{k]\fj^i)\^ Note that since D and V commute, we can simultaneously impose equations 
(I6.23P and (16.241) . In this case, we write f{p) for a source of momentum p. 

We now construct the form of these solutions in the fiat space limit. To do this, we briefly 

At finite N, this condition has various correction terms. Indeed, whereas the 50(5) generator 
V is Hermitian, P is nilpotent. In the flat space limit, however, this is not much of an issue. 
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review the construction of solutions to the free wave equation in commutative twistor space. 
Solutions to the helicity h free field wave equation on 4D spacetime are conveniently specified 
by the Penrose transform. The basic idea is to look for elements of H^{FT' ,0{2h — 2)), 
where PT' = CP'^ — CP^ is projective twistor space with the line at infinity deleted. Given a 
cohomology representative /p(aj, ir) on PT', we can via the twistor equation u = ixn obtain 
a representative on correspondence space x CP^, with coordinates {x"'°',Sa)- Observe 
that this is not an arbitrary section on the correspondence space; It satisfies the condition: 

X''daaU{x,s) = (6.25) 

The Penrose transform amounts to a contour integral over the CP^ factor, resulting in a 
4D space-time field: 

0(x) = (f{sds)lpix,s) (6.26) 



where here, we have specialized to the case where fp is a degree zero (0, l)-form. Other 
helicities are covered by including factors of tt or d/du acting on the integrand fp. The 
plane wave solution is: 

fp{x, s) = exp(ip ■ x)5{{s\)) (6.27) 

where the delta function is a (0, 1) form of specified homogeneity. The Fourier transform 
of equation fl6.26p then provides a helicity h momentum eigenstate. 

Let us now turn to the fuzzy setting. The analogue of the delta function 6{{sX)) on a 
fuzzy CP^ is the operator |A)(0|. Summing over a basis of position states, much as in our 
discussion of the projection to Minkowski space, the corresponding operator f{p) is: 

f{p) = J d'^x exp{ip ■ x)\x, X){x,0\ (6.28) 

where Pa = A^Aq. Observe that this operator formally satisfies the operator equation 
7r"[-Pda, /(p)] = 0, which is the analogue of equation f l6.25p . To obtain the corresponding 
plane wave operator V{p), we need to integrate equation fl6.23p . Here we can make use of 
the CP^ propagator fl6.16p . which satisfies f l6.17p . We thus arrive at the following definition 
of the plane wave operators: 

V{p) = j d^x e'P-''\x,\){x,\\. (6.29) 

The integral expressions f l6.28p and f l6.29p should both be viewed within the context of the 
scaling limit (16. ip . This means that the integral runs over a local patch near x = 0, while 
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the momentum p is scaled accordingly (see figure [T]). 

The generalization to the supersymmetric case is straightforward. The plane wave states 
then also depend on anti-commuting variables ( 

V{pX) = y"rf^l^a;e'P"+'^-<^^>|x,^,A)(x,^,A| (6.30) 

The components of the = 4 supermultiplet V provide states of different helicity in the 
4D theory. For the A/" = 4 vector multiplet, the bottom component constitutes a plus 
helicity gluon and the top component is a minus helicity gluon. For A/" = 4 gravitons, the 
plus helicity and minus helicity gravitons sit in two different supermultiplets, with the plus 
helicity mode at the bottom of its multiplet, and the minus helicity mode at the top of its 
multiplet. As noted in [H], MHV scattering amplitudes in N' = A superspace lead to a 
kinematic factor of (nl)^ for minus helicity states of momenta pi = AiAi and p„ = A„A„. 
To avoid clutter, we shall often leave implicit the integration over A/" = 4 superspace. 



6.4 MHV Gluon Scattering 

We are now ready to compute amplitudes. We first consider the MHV gluon scattering 
amplitudes [26]. The computation is similar to that in [H] and to the twistor string theory 
calculation in [12]. As discussed in section 3.2, the gaussian matrix model comes with a 
natural set of currents J^iV), given in fl3.17p . that describe the coupling to a gauge field A, 
defined on non-commutative CP"^. Here the V are linear maps acting on the CP^ Hilbert 
space, and should be thought of as the asymptotic wave functions of the gluon states. A 
supermultiplet of a gluon state with color charge and supermomentum (p, () is described 
by 

J^{p,0 = Ti (r^V{pX)QDQ'^ (6.31) 

where V{p, () is the momentum eigenstate defined in (16.301) . 

The scattering amplitude is now directly given by the matrix model expectation value 

iMl,...,n= (jl{pi,Cl)-Jn{Pn,Cn)) • (6.32) 
\ / MM 

As explained in section 3.2, performing the matrix integral is trivial: the vertex operators 
( I6.3ip represent the response to a simple field redefinition SyQ = VQ, 6vQ = 0. Performing 
resulting Wick contractions immediately leads to the following expression for the color- 
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stripped subamplitude 

= Tr (V{p,, Ci) ■ ■ ■ V{pn, Cn)) (6.33) 

Here the trace is over the Hilbert space CP'^'^ of the non- commutative supersymmetric 
Hilbert space. 

The rest of the calculation is equally straightforward. Inserting the definition f l6.3Up for 
the plane wave operators, we first evaluate the amplitude in position space 

n n r4|8/ \ 

1 1 

1=1 1=1 * ' 

where i = is identified with i = n. Here we used the result fl6.18p for the overlap. Fourier 
transforming back to momentum space, we obtain 

-^1&)m^ (''•^«) 

which we recognize as the color-stripped contribution to the gluon MHV amplitude. Let 
us note that the form of this amplitude is fixed up to an overall multiplicative constant. In 
particular, there is an overall finite factor of 5(0), which is connected with a short distance 
cutoff. Here, this simply specifies a reference energy scale for the external momenta pi of 
the amplitude. 



7 MHV Graviton Scattering 

In this section we show how MHV graviton amplitudes are computed in the matrix model. 
MHV graviton scattering amplitudes involves a pair of minus helicity gravitons (labeled by 
1 and n) and an arbitrary number of ingoing plus helicity gravitons (numbered 2 to n — 1). 
In this section we demonstrate that this amplitude is given by the expectation value of 
gravitational currents in the matrix model: 

iMmnW = (fl{pi,Cl)T+{p2) ■ ■ ■T+{Pn~l)fn{PnXn)) (7.1) 
\ / MM 

The twistor realization of this calculation builds on Penrose's non-linear graviton construc- 
tion of ( ant i-) self dual space time backgrounds. Collectively, the plus helicity gravitons can 
be thought of as representing a selfdual background geometry, on which the minus helicity 
graviton propagates. Upon reversing the momentum of one of the minus helicity gravitons 
(say, the one labeled by n), the MHV scattering amplitude represents [33] the process of 
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Figure 2: The MHV graviton scattering amplitude represents the process of a minus helicity 
graviton reflecting off a self-dual background, built up by the metric fluctuations with plus 
helicity. The scattering produces a helicity flip. In the matrix model this amplitude is 
computed in a small "Minkowski patch" . Figure modifled from [33] • 

an incoming minus helicity graviton bouncing off the selfdual background geometry, while 
undergoing an helicity-flip. The MHV amplitudes arise by expanding out the background 
fleld in terms of linearized perturbations around flat space-time. 

Graviton excitations are associated with deformations of space-time. We have already 
identifled a natural class of glik^) generators in section 3.2 given by equation fl6.22p : 



in the obvious change of notation. Here v-V = v"'°'Paa + VaaK°''^ is the (complexifled) 5*0(5) 
generator, that in the local Minkowski patch around a; = acts like a translation generator 
in the direction v"''^. These are momentum eigenstates connected with translations. Note, 
however, that it also contains a component proportional to the conformal boost generator. 

Before explaining the geometric content of these currents, let us flrst specify the form 
of the tensor v'^"' appearing in Ty. Here we will slightly simplify our task. Since we are 
anticipating an interpretation in terms of graviton states, we can carry over the standard 
space-time treatment of the associated polarization tensors. Of the four possible polar- 
izations w"", we can eliminate two due to spin 2 gauge invariance: linearized space-time 
diffeomorphisms shift f "'^ with an amount proportional to the momentum p"''^. Such longi- 
tudinal modes decouple from the amplitude (we will make this explicit shortly). This then 
allows us to choose the transversality condition VaaP"""" = 0. The remaining two polarizations 



T,{p,0 = {v-V)V{p,C) 



(7.2) 
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split up into a positive and minus helicity component via 



(/iA) = 1 



(7.3) 



aa 



[/iA] = l 



(7.4) 



This choice of the dual two component spinors Ha and fia is not unique, due to the residual 
on-shell gauge invariance /ia — /ia + Aa and fia f^a + Ad- We will make a specific choice 
for fia and fia later on. Note that the transversality condition p°'°'Vaa = has the natural 
consequence that, in the fiat space limit, the generator v"-"'Vaa v°'°'Paa commutes with 
V{p, (). To denote the two types of helicities, we will continue to write T+{p) or just T{p) 
for the plus helicity gravitons, while we will write T-{p,() or just T{pX) for the minus 
hehcity gravitons. 

It is instructive to look at the leading order form of the plus helicity graviton currents 
in the local Minkowski region near u = 0. In this fiat space limit, the polarization operator 
of the plus helicity gravitons takes the form 



where we used that 'n'"'V{pk) = X1^V{pk) and (/ifcAfc) = 1. Here the invariance under shifts 
/^a — ^ /^a + Aa is made manifest. Due to the non-commutativity, equation f l7.5p represents 
a holomorphic derivative along the Ua directions. The vertex operators of the plus helicity 
gravitons can thus be viewed as living in the holomorphic tangent space to twistor space. 
Moreover, we learn that the insertion of the plus helicity generator T^{p) represent small 
transverse deformations of the local twistor line, via an infinitesimal diffeomorphism that 
shifts Ua by an amount proportional to XaV{p). Similar considerations hold for the minus 
helicity generator T_{pX)- 

Having specified a class of gl{k]si) elements, we next specify the corresponding vertex 
operators. As we have already mentioned, the gl{ki^) x gl{kN) gauge symmetry are associ- 
ated with the bulk gauge field A and the compensator A. As opposed to the case of MHV 
gluon scattering, in MHV graviton scattering the two helicity modes fill out distinct A/" = 4 
supermultiplets. This means that we should anticipate that the presentation of their vertex 
operators can be different. In an MHV amplitude, we expect the plus helicity modes to 
deform the local geometry, while the two minus helicity modes specify asymptotic data. 
This motivates our specification of the vertex operators: 



{v^rPaa=M[~XkUJ^ = [~Xuj^ 



(7.5) 




(7.6) 
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Here, the plus helicity currents 7+ act via the algebra of vector fields, and embody the 
result of performing an adjoint gauge variation 

(5g= [T+(p),g] ; 5g = o (7.7) 

Geometrically, these transformations build up a selfdual background off of which a minus 
helicity graviton can scatter. The minus helicity gravitons specify flow into and out of the 
patch near the south pole, and transform in a dual representation on which the vector flelds 
can act. 

This identiflcation makes direct contact with the continuum twistor description of gravi- 
ton scattering amplitudes [33]. In Penrose's non-linear graviton construction [29], a selfdual 
background space-time is converted into a twistor space with a deformed complex structure. 
The deformed space is still described by the same coordinates as usual twistor space, and 
TVa still represents a holomorphic coordinate along the twistor lines. The coordinates u"" are 
however no longer holomorphic in the distorted complex structure. 

To bring out this geometric picture, let us introduce the self-dual gravity background 
given by the formal exponentiation of the plus helicity metric fluctuations 

n-l 

exp {T+{h)) = H exp {hT+{pk)) (7.8) 

k=2 

where is the inflnitesimal amplitude of the k-th mode. Each individual graviton amounts 
to a small chiral GL{ki^) gauge transformation, which acts via commutation as in fl7.7p . 
This transformation represents the non-linear graviton background. We can formally write 
the action of the gaussian model propagating in a self-dual space-time background as 

^,(g,g) = Tr(e^+ wg^g) (7.9) 

where the expansion of the exponential acts by all possible orderings of nested commutators 
on g. As explained, the MHV graviton scattering amplitude can be obtained by considering 
the propagation of a minus helicity graviton on top of this self-dual background. 

Consider next the vertex operators for the minus helicity gravitons. As opposed to the 
case of MHV gluon scattering, in MHV graviton scattering the two helicity modes flU out 
distinct A/" = 4 supermultiplets. This means that we should anticipate that the presentation 
of their vertex operators may be different. As we have seen, the plus helicity modes are 
identifled with holomorphic vector flelds that deform the local geometry. Hence, following 
the continuum twistor description, it is natural to view the minus helicity modes as being 
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in the cotangent space, which are transformed by this algebra of vector fields. We will 
identify the minus helicity graviton with the following currents 

f^{psXs) = TT(f_{psXs)QDQy s = l,n (7.10) 

with T_ as in equation (17. 2p . with polarization tensor appropriate for a minus helicity 
graviton. The momentum associated with the minus helicity gravitons flows in and out of 
the local Minkowski patch, as indicated in fig. [2J 

The geometric interpretation of (17.101) as associated with cotangent factors arises as 
follows. The polarization factor of the minus helicity graviton, say graviton 1, contains the 
contraction of (f f)aa with the conformal boost generator K°''^, 

(v^UK'"' = {X^n^)[^^^cu]. (7.11) 

When we compute the amplitude, we will see that this operator naturally pairs with the 
holomorphic tangent vectors (17. 5p . in the precise way that one would expect from an element 
of the cotangent space. 

7.1 Evaluation of the Graviton Amplitude 

We will now calculate the matrix model expectation value (17. ip and show that it reproduces 
the MHV amplitude (II. 4p . Our notation for the plus and minus helicity currents is Vk = 
{v^ -V) , for k = 2, ..,n — 1, Vs = {v~ -V) for s = 1, ra, and Vi = V{pe) for all i = I, .., n. 

The calculation is remarkably straightforward. The basic form of the combinatorics is 
similar to the computation in [33] and also has some overlap with the proposal of |l2]. Using 
the exponentiated form of the plus helicity gravitons as representing a self-dual background, 
performing the Wick contractions produces the following amplitude 

Tr(fi-exp(r+(/i)) -f^) (7.12) 

This expression exhibits the geometric interpretation of the amplitude as the propagation 
of graviton 1 through the self-dual background, emerging as an outgoing graviton n. See 
figure [3] for a depiction of how the non-linear graviton deforms the local twistor geometry. 

One could in principle try to do the calculation with the fully exponentiated background. 
Instead, let us expand the self-dual background into the individual graviton modes. Ex- 
panding out each contribution from T+, we obtain a set of nested commutators for our 
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Figure 3: Depiction of the non-linear graviton (red arrows) in fuzzy twistor space. The 
dashed hne indicates a CP^ with respect to a fixed complex structure. The curved line 
indicates the deformation due to the non-linear graviton. Each red arrow corresponds to 
a "bit" of the background and the cumulative effect generates a finite complex structure 
deformation. See [33] for a related depiction of the non-linear graviton. 

subamplitude: 



The continuum theory interpretation of this form of the amplitude is hopefully clear; The 
nested set of commutators correspond to the action of Lie derivatives which act to the left 
on Ti. Later, we will see that Ti can be viewed as transforming in the cotangent bundle. 
The successive Lie derivatives then return a modified cotangent vector which dots with the 
(dualized) cotangent vector T„. The full amplitude is then given by summing over P(2,...,n-i), 
that is, all permutations of the plus helicity gravitons0 Finally, in the above expression, 
we have suppressed an overall prefactor of k = ^/16nGN for each graviton vertex operator. 
The precise value of Newton's constant cannot, however, be fixed by just MHV graviton 
scattering. Instead it requires a discussion of how the gaussian matrix model embeds inside 
a more complete framework of the type discussed in the companion paper [16], that include 
modes that propagate between different twistor lines. 

We would like to evaluate the subamplitude while taking the scaling limit (16. ip . 

As opposed to the case of gluon scattering, implementing this scaling limit for graviton 
scattering is slightly more subtle, because the currents are themselves related to infinitesimal 
translations. Indeed, the correlator we are considering is technically defined with respect 
to an S"^ geometry, while the scattering amplitude is defined by taking a scaling limit in a 

^^Returning to equation (j7.1|) . we should note here that Wick's theorem also produces additional 
orderings of the currents where Ti and T„, are not adjacent to one another. These do not correspond 
to graviton scattering, but rather describe a graviton which would be absorbed into the self-dual 
background. In the flat space limit, all of these contributions vanish simply due to conservation 
of helicity and momentum. 




(7.13) 
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small patch near the south pole. Roughly speaking, we need to take into account the fact 
that whereas in Minkowski space we can have a momentum flow in from inflnity, on the 
compact S*^ geometry, this requires introducing a source and a sink. 

A convenient way to capture the large radius behavior of the correlator is as follows. The 
flat space limit is obtained via insertions of the operator 1m introduced in equation fl6.6p . 
For the most part, these insertions can be ignored. However, the insertions involving the 
minus helicity states must be treated with more care, because they can move the domain of 
definition which we are zooming in on. Focussing on the insertion between the two minus 
helicity graviton currents, we can insert a factor of 1m into the subamplitude: 



which implements the flat space limit. Recall from equation (16. 6p . the projection to the 
patch is obtained by a sum over \x, p){x, p\. Here we adopt angular brackets to indicate we 
are now working about a patch of Minkowski space. To leading order in the flat space limit, 
the translation generators annihilate this state via right multiplication: \x, p){x, p\Paa = 0. 

We now consider the subamplitude ... ,i in this scaling limit. To simplify the expres- 
sion, we flrst expand out each successive commutator. Doing so, we see that for the plus 
helicity generators Tj, \x,p){x,p\Ti = 0. In other words, each nested commutator can be 
replaced by ordinary matrix multiplication. To represent the flow of the momentum into 
the patch, we keep Vi to the left of Vi and we take Vn to the right of Vn, which can be 
done since [Pn, Vn] — in the flat space limit. Our new task is therefore to evaluate: 



where all operators now compose via ordinary matrix multiplication. In this expression we 
have made the substitution Tj — )■ PiVi for all the plus helicity generators, as appropriate 
in the flat space limit. This form of the subamplitude has the clear interpretation of the 
minus helicity graviton flowing into the local patch near the south pole, scattering off a 
self-dual background, and then flowing out of the local Minkowski patch. 

The basic plan of the calculation is therefore very simple: we flrst transport the k = 
2, ...,n — 2 plus helicity factors of Pk all the way to the left. Along the way we collect all 
contributions of the commutators with the plane wave operators V{pi) 




(7.14) 




(7.15) 




(7.16) 
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Note that in the hmit we are considering, these momentum factors give the dominant con- 
tribution: the contributions from (x,p|Pfc and from moving past Vi are both neghgible. 
So the leading order result is: 

n-2 

{x,p\ViV^^2...Vn-2Pn-lVn-lVnVn\x,p) X{vt ■ [pi + . .. + Pk-l) ■ (7.17) 

k=2 

Note that alternatively, we could have allowed all of the P's to pass to the right. This 
would result in factors of v'^ ■ {pk+i + ... + Pn)- Via conservation of momentum, this is the 
same as the factor obtained in equation f l7.17p . In [33] a simplifying gauge choice was taken, 
where fik A.„ for all k = 2, ...,n — 1. In this gauge, if we pass Pn~i to the left so that 
it sits just to the right of Vi, the leading order contribution from the commutators fl7.16p 
vanishes because, via momentum conservation, this term is proportional to v^_i-pn = 0. 

So to evaluate the Pn-i contribution, we first commute the oj^ subfactor to the left, 
where it hits the uja oscillator in the Ki term in Vi = Pi + Kii~^. Using [w^, Ua] = e^h, and 
that at large A^, we can replace 'K'^°"Kb = N6°'b, we obtain for the subamplitude 

N ""^ 

- {x,p\V,...VnVn\x,p) K--^,tl) llv+-{p, + ...+Pk-l). (7.18) 

k=2 

Notice the presence of the additional prefactor N/i"^ ^ V^pi- The amplitude thus remains 
finite in the scaling limit fl6.ip . Including the overall normalization factor k", the amplitude 
then scales expected for the MHV graviton amplitude. 

The last step is to take into account the flow of momentum out of the patch, as indicated 
by the factor of P„. The leading order behavior is obtained by taking P„ — P„, which 
induces a small translation on the ket state \x,p) — )■ |a; + fn,p). This has the consequence 
that it shifts the domain of integration in equation (16. 6p . The effect of the shift is given 
by the value of the subamplitude with \x,p) held fixed but with P„ acting to the left on 
Vi. This amounts to the substitution Vi — (?;„ ■ Pi)Vi, and reflects the fact that the minus 
helicity graviton backreacts on the patch as it flows into and then out of the patch. 

Now that we have taken into account the gravitational "color factor", the remaining 
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manipulations are the same as for gluon scattering: 



n-2 



^,...,n = «"-' Tr(ri...K) {vi-vU){v--Pi) l[v+-{p^ + ...+pk.,) (7.19) 

k=2 

n-2 c-4|4 (^1 ■ ^'n-l) ('''n ' Pi) TT + / . . x r,n\ 

= « 5' X TTT^^T ii^fc -(Pl + ••• +Pfc-l) (7.20) 

^ ^ k=2 

where the factor of 5^^^ enforces conservation of momentum on the TV = 4 superspace. Here, 
C{n) = (12)(23) ■ ■ ■ (n — ln)(nl) denotes the Parke- Taylor denominator. 

To complete the computation, we plug in the values of the reference twistors used in 
our evaluation of the current algebra. For the plus helicity polarization tensors, we have: 

= 1^ (7^21) 

for k = 2, ...,n — 1. The dual spinor (17.211) automatically satisfies (fik^k) = 1- Note that 
this choice implies that (fikn) = 0. The minus helicity polarization tensors are: 

Inserting this explicit form of the polarization tensors, we obtain: 

V^^V [/3n](n-ln)(ln)' C(n) Al (kn) ^ ' 

with C{n) = (12)(23) ■ ■ ■ (n — ln)(nl) the usual Parke- Taylor denominator. Rather re- 
markably, this is identical to the form of the subamplitude found in [33]! At this point, 
the combinatorics and manipulations to show that this expression reproduces the original 
BGK result (11.41) are the same as in [33]. The full amplitude is now obtained by summing 
over all permutations of the plus helicity gravitons 2, ...,n — 1. As explained in [33], note 
that although the amplitude appears to depend on the gauge choice /3, once we sum over all 
possible permutations, this dependence drops out, as it must. For formulae with a manifest 
Sn-2 permutation symmetry of all plus helicity gravitons see [13]. 
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8 Conclusions 



In this paper we have studied the properties of a gaussian matrix model formulated on 
fuzzy twistor space. The model possesses a number of symmetries which at large N have 
a direct interpretation on a 4D space-time. Identifying a natural class of currents for the 
theory, we have shown that correlators of the matrix model reproduce both MHV gluon and 
graviton scattering amplitudes. This provides non-trivial evidence that the twistor matrix 
model of [121 HT] correctly describes 4D physics, and moreover, contains a gravitational 
subsector [TB] . 

Our discussion has been limited to MHV amplitudes. To recover more intricate ampli- 
tudes, we expect that the details of bulk physics in twistor space will be important. This 
is in keeping with the qualitative picture in twistor string theory that MHV gauge theory 
amplitudes localize on lines of twistor space [12]. Let us note that the gravity sector adds 
an interesting twist to this; These amplitudes essentially localize on a complex line, but one 
which has been deformed by the background gravitons. To make further contact with 4D 
physics, it would be interesting to study more general types of amplitudes in the framework 
proposed in [T6] . 

Since we have an explicit finite matrix model, we can also study various cor- 
rections to the continuum theory result. In the flat space limit, the corrections of physical 
relevance are controlled by the small parameter ipi. More broadly, one can also consider 
correlators on a finite size S"^ which may provide a concrete means to compute (via analytic 
continuation) in-in correlators on de Sitter space. 
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A Twistors and SO{5) 

In this Appendix we review some aspects of the infinity bitwistor for an S^. See [TB] for 
a somewhat expanded discussion. With the help of the infinity twistor, we can define 
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normalized six component space-time coordinates x"^, and their duals Xa/3, via 

These satisfy the relations x°'^ x^'^ = x"^, Xa'yX'^ p = Xap,- The dual coordinates Xap act as 
projection matrices onto the twistor line associated to the space-time point X: the twistor 
line equation 

x^fiU^ = (A.2) 

is solved by all points U = aZ + bW on the twistor line through Z and W. 

The dual space-time coordinates x"^ can be parametrized as with the help of five coor- 
dinates constrained to live on an of unit radius, via 

Xap= \ , ^ ^ \, y Va = 1 • (A.3) 

V i^'^-y^KbJ 

Here = \{y^o''^)b with = (a*, — ^1) the usual Pauli matrices. With this parametriza- 
tion, the four component twistor line equation ( \A.2\i reduces to the standard two-component 
twistor line equation with 

The fiat space limit amounts to zooming in on the south pole region of the near ~ 
— 1. In this limit, the remaining four coordinates become identified with fiat space 
coordinates x'^. 



B Fuzzy Cauchy 

In this Appendix, we give an alternative construction of the dual position state (A| in 
'H*pi(A^) that satisfies the inhomogeneous coherent state property (15.191) . 

The dual Hilbert space H^pi (A^) is + 1 dimensional. So let us introduce a convenient 
position eigenbasis, by picking a preferred set of A^ + 1 positions Zp, with p = 0, .., N, for 
which we will then construct dual position eigenstates. We will then later define general 
position eigenstates for other locations via a suitable version of analytic continuation. A 
natural choice is to take the special positions Zp to lie on the unit circle 

Zp'^^ = 1, so that Zp = exp(^ ^"^^ (B.l) 
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We now introduce the discrete set of candidate dual position eigenstates, via 



TV 



n=0 



These form a complete basis of 'H*pi(A^). Moreover, they are indeed almost dual position 
eigenstates. The violation of the coherent state condition is localized at the two poles: 

{zp\ (tti - Zp'K2) = (0| - (oo| (B.3) 

where (0| and (oo| = (A^| denote the dual north and south pole state in H^^i{N + 1). 

To extend the definition of the dual position states to other points away from the unit 
circle, we will use a discrete version of Cauchy's formula. We introduce the following 
notation 

1 ^ 

p=0 

Indeed, we claim that the continuum limit of the right-hand side amounts to performing 
a standard complex contour integration along the unit circle. As a first trivial check, we 
compute the basic Cauchy integral 

^ ^.(0 (B.5) 



N+l 



In the large limit, the function on the right hand side reduces to a step function: it 
is equal to 1 on the northern hemisphere where |^| < 1, and vanishes on the southern 
hemisphere where |.^| > 1. So at large N we can replace 6*^(0 by 

1 for 1^1 < 1 

0(0 = { (B.6) 

for 1^1 > 1 

This is indeed what one expects from a discretized residue theorem: when ^ traverses the 
equator, it sneaks between the holes in the discrete contour and escapes. We will return 
to this special feature of the discretized contour integration momentarily. 



^^On the southern patch, when |^| < 1, we need to use the dual coordinates A = ^ ^ and 
p = . The Cauchy formula then becomes §^ ^ \^ = 6^ (A) where §^ = §^ 
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With this new tool in hand, we now define dual position states for any position ^ via 

(el = / i^p\ (B.7) 



This provides a well defined state for any value of ^, except for solutions to = 1. The 
state (^1 almost satisfies the dual coherent state requirement. Upon inserting the definition 
flB.7l) . while using equations fIB.Sp . (IB.SP and the formula {zp\7r2 = (oo|, one finds after 



a straightforward calculation 

(el(vri-evr2) = ^..(0(01 +^.(r')(oo| (B.8) 

with (0| and (oo| the dual south and north pole state in 'H'^^i^N + 1). In the large limit, 
equation fl5.19p reduces to 

f (0| for|e|<l 
m^i~^^2)={ (B.9) 
[ (oo| for 1^1 > 1 

This is our desired intermediate result. It shows that the dual coherent state (.^| are position 
eigenstates, modulo a source term localized at the corresponding pole. 

Finally, let us compute the overlap between the position eigenstates. A simple calcula- 
tion shows that the dual basis {zp\ pairs with the position states ( l5.1Up via 

(^,|A) = -J— (B.IO) 

i — Zp A 

Combining this result with (lB.7p . we obtain that 



This equation reveals, as expected, that ^ and A are reciprocal affine coordinates. Note, 
however, that the pole in the denominator is spurious. The step functions do not allow ^ 
and A to be located on the same hemisphere: whenever they do, the numerator vanishes. 

The standard way to overcome this obstacle is via analytic continuation. Consider the 
overlap (^|A) and let Ai be the reciprocal coordinate to .^i = Aj"^. We wish to define the 
south patch state (Ai| via analytic continuation of the north patch state However, 
here we meet a subtlety. At infinite A^, the step functions are non-analytic at the equator, 
reflecting the jump over the Cauchy contour. At finite A^, on the other hand, the step 
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X 



•2 

X 



Figure 4: Evaluation of correlators of the chiral boson on a fuzzy CP^ is achieved by working 
in terms of a basis of states which vary holomorphically as a function of the complex plane 
parameter A. In the figure, two patches of the CP^ are indicated by the dashed circular 
contour. Correlators in different patches of the CP^ are computed by first evaluating in one 
patch in the large N limit, and a subsequent analytic continuation. 



functions 6'jv(0 perfectly analytic, reflecting that the contour is just a discrete set of 
points. The normal strategy of analytic continuation would amount to pushing the location 
of the contour. Hence we can proceed in two ultimately equivalent ways. We can either 
first take the large N limit, and then analytically continue. Or we can define analytic 
continuation by suitably deforming the discrete contour. In practice this means that the 
location Ai = A2, where we expect the pole to occur, must lie somewhere on the discrete 
contour. At very large A^, this approaches the continuum prescription, except very close to 
the pole. At infinite A^, the two prescriptions coincide. 

The upshot is this: we define the state (Ai| such that its overlap with IA2) is given by 
the analytic continuation of (IB. lip , starting from the region 6{^) — 6{X~^) = 1. Hence in 
the strict large A^ limit, we have 

This completes our construction of the state (A|. 
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